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The adsorption geometry and electronic properties of a zinc-phthalocyanine molecule on a Cu(lll) 
substrate are studied by density functional theory. In agreement with experiment, we find remarkable 
distortions of the molecule, mainly as the central Zn atom tends towards the substrate to minimize the 
Zn-Cu distance. As a consequence, the Zn-N chemical bonding and energy levels of the molecule are 
significantly modified. However, charge transfer induces metallic states on the molecule and therefore is 
more important for the ZnPc/Cu(lll) system than the structural distortions. 

Understanding and modifying materials properties for technological applications has always been a prime 
interest of materials science. To develop new functional materials for nanoscale electronics, it is important 
to accurately characterize and manipulate the quantum state of molecules adsorbed on surfaces. Organic 
semiconductor devices have attracted much attention due to various potential applications^ "I Metal phthalo- 
cyanine (MPc) molecules are important for organic electronic devices because of their good chemical and thermal 
stability, less complicated synthesis, and tunable physical properties by substitutional flexibility^'^. They provide 
relatively stable chemical environments for the central metal atom^. In addition, the molecules can be employed as 
building blocks for a wide range of systems such as organic thin film transistors, field effect transistors, gas 
sensors, photovoltaic cells, and organic light emitting diodes^. 

It is observed that the device functionalities of organic molecules can be significantly improved by controlled 
selective adsorption on different substrates Of course, the electronic properties of the molecules are modified 
by interaction with the substrate. In particular, charge transfer between substrate and molecule modifies the 
molecular band gap^'^^'^l Recently, Zhang and coworkers have investigated the adsorption of MPc molecules 
(M = Fe, Co, Ni, Cu, and Zn) on Au(lll) and have established the most stable configurations^. They report on 
significant differences in the electronic properties of MnPc, FePc, and CoPc molecules. Baran and coworkers have 
presented a combined experimental and theoretical study of adsorption of CoPc, PbPc, and SnPc on Ag(lll)^^. 
They find a strong interaction of the central Pb atom with the Ag surface, leading to a reduction from Pb^^ to Pb°. 
Among the MPc molecules, ZnC32N8Hi6 (ZnPc) has become a kind of model system for fundamental studies. It 
shows most of the basic features of the Pc molecule, because the influence of Zn on the electronic structure is small. 
Numerous optical and electrical applications of ZnPc in organic electronics have been discussed in the literature^'^^. 

The supposition that aromatic molecules without electronic coupling are physisorbed on metal surfaces is not 
acceptable for all systems. Though there is no electronic coupling and, hence, no covalent bonding in the 
pentacene/Cu(100) system, it has been observed that the electronic structure of pentacence is modified signifi- 
cantly. This modification is due to the formation of mixed molecule-metal states^^. Formation of interface states 
and dipoles is a central issue for active organometallic interfaces They determine the energy level alignment 
near the Fermi energy (Ep) and thus the performance of the device. Key factors in the creation of interface states 
and dipoles are the distance between molecule and substrate and the transfer of charge between them^^'^°. The 
adsortion of a molecule can modify the work function of the substrate and thus the interface dipole, which 
determines the carrier injection barrier at organometallic interfaces^^'^^. 

Yamane and coworkers have used a combination of the x-ray standing wave technique and angle-resolved 
photoemission spectroscopy for ZnPc on Cu(lll) to analyze how the adsorbate configuration affects the chem- 
ical bonding at the interface and its electronic properties^ \ In addition, Gerlach and coworkers have shown 
experimentally that the adsorption distance on a Cu substrate is small because of the relatively low work function 
(as compared to Ag and Au, for example), pointing to a strong interaction^^. As a consequence, charge can easily 
be transferred between the substrate and adsorbed molecule to modify the electronic properties. Despite many 
experimental results on the adsorption of organic molecules on metal surfaces, an adequate theoretical descrip- 
tion remains challenging^^. 
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Figure 1 | Free standing ZnPc molecule. 

In this article, we determine the adsorption geometry and elec- 
tronic structure of an isolated ZnPc molecule on Cu(lll). Going 
beyond the experiment, we are able to describe and explain the con- 
sequences of the interaction with the substrate in detail. On the one 
hand, we demonstrate that the Zn-N chemical bonding within the 
molecule is modified such that the energy levels are altered. On the 
other hand, a strong transfer of charge within the molecule leads to a 
metallic nature. This fact indicates that the charge transfer governs 
the properties of this organometallic interface rather than structure 
effects as proposed in Ref 21. 

Results 

In order to investigate the adsorption behavior and electronic prop- 
erties of ZnPc on Cu( 111 ), it is necessary to first study a free standing 
ZnPc molecule. Figure 1 shows the schematic structure of the mole- 
cule, which is composed of a central Zn atom and four isoindole 
groups bridged by four additional N atoms, having planar D4h sym- 
metry^^. To analyze the electronic structure of the free standing ZnPc 
molecule, the total and partial DOSs are calculated and presented in 
Fig. 2. It is clear from the total DOS in panel (a) on the left hand side 



of Fig. 2, which is obtained within the GGA, that the highest occupied 
molecular orbital (HOMO) and lowest unoccupied molecular orbital 
(LUMO) are due to the Pc skeleton, with small contributions of Zn to 
the HOMO. These two highly delocalized orbitals are responsible for 
most of the physical properties. The state near Ep is due to the C 2py 
orbitals, while the states around — 2 eV are a combination of C Ip^, 
2px and N Ip^, Ip^ orbitals. Similarly, the LUMO is due to both C Ipz 
and N Ipy orbitals. The group of peaks around —4 eV is mainly due 
to the C 2pz orbital, with small contributions of the N 2py orbital. In 
addition, the states in the energy range from — 8 eV to — 5.5 eV trace 
back to the C 2p^, 2px and Zn 3d orbitals, with a very small contri- 
bution of the H Is orbital, see panels (b), (c), and (d) in Fig. 2. We 
observe />(i-hybridization between the C 2p and Zn ?>d orbitals in this 
energy range, where Zn contributes more at low energy. This fact 
indicates that the presence of Zn at the center of the molecule does 
not influence the basic electronic structure of the Pc skeleton, reflect- 
ing the suitability of ZnPc as model system for fundamental studies. 

We have also calculated total and partial DOSs of a free standing 
ZnPc molecule using the hybrid functional approach (PBEO), see the 
right hand side of Fig. 2. The 1 eV peak in the LUMO (C 2p^ and N 
2py States) is shifted to around 2 eV, while the small Zn 'id peak below 
Ep shifts to lower energy in the case of the hybrid approach. Thepii- 
hybridized and H states reveal a similar behavior. We note that our 
GGA DOS of free standing ZnPc shows significant improvements 
over previously reported results in Ref 28, where the local density 
approximation has been employed, and seconds Ref 5. The PBEO 
approach yields a band gap of 1.93 eV, in excellent agreement with 
the experimental value^^. In addition, the photoemission spectra 
reported in Ref 30 are well reproduced, particularly the small feature 
at -0.24 eV which appears both in our GGA and PBEO data. While 
the states found in the experiment around —2.5 eV could not be 
reproduced in previous GGA and GW calculations, they clearly 
reappear in our case. This difference may be related to the C pseu- 
dopotential, as these are C states. 

Before discussing the electronic structure of the ZnPc/Cu(lll) 
system, we first address the adsorption geometry (relaxation effects). 
To this aim, we investigate three adsorption sites: top, bridge, and 
hexagonal closed packed (hep). We also consider the angle between 
the wings of the ZnPc molecule and the substrate. Therefore, we deal 
with six configurations in total: top, top angle, hep, bridge-I (B-I), 
bridge-II (B-II), and bridge-Ill (B-III). Figure 3 gives top views of 
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Figure 2 | Total and partial DOSs (per atom) of a free standing ZnPc molecule obtained by the (left) GGA and (right) PBEO approaches. 
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Figure 3 | Adsorption configurations of ZnPc on Cu(lll). (a) top, 
(b) top angle, (c) hep, (d) B-I, (e) B-II, and (f) B-IIL 

these six configurations, which we relax with and without taking into 
account the vdW interaction to obtain the adsorption energy {Bads)- 
The results are given in Table 1. We find the top configuration to be 
most stable, in contrast to Ref 5 where the top angle configuration 
was found to be favorable for the ZnPc/Au(l 1 1) system. This differ- 
ence reflects the different substrates. Higher adsorption energies for 
other configurations do not exclude metastable states at elevated 
temperatures, in particular the hep configuration is close in energy. 

After relaxing the system, we obtain a significant distortion of the 
ZnPc molecule and find substantial differences in the distances 
between molecule and substrate for calculations with and without 
vdW interaction. The average distances with vdW interaction are 
3.22 A, 3.09 A, 2.98 A, and 2.67 A for H, C, N, and Zn, respectively. 
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Without vdW interaction we obtain 3.12 A, 2.98 A, 2.77 A, and 
2.42 A, respectively. Therefore, in general, the distances are larger 
under inclusion of the vdW interaction, showing that the vdW 
method does not yield good results for the present system. An over- 
estimation of the distance between an aromatic molecule and Cu 
substrate with vdW interaction has also been reported in Ref 31 
and agrees with the conjecture of Ref 5 that the vdW method usually 
has drawbacks for molecule-metal interfaces. Our calculated dis- 
tances obtained without vdW interaction deviate from the experi- 
mental values by 19.6% for C, 8.6% for N, and 6.6% for Zn. RuraH and 
coworkers have reported similar deviations for a perylene derivate on 
Ag(lll)^^. Toyoda and coworkers have argued that the overestima- 
tion is due to the more repulsive exchange part of the generalized 
gradient approximation^^ 

An interesting outcome of our calculations are the significant 
differences in the distances of the different atoms of the molecule 
from the substrate. They reflect a strong distortion of the ZnPc 
molecule (non-planar configuration), where the Zn atom is pulled 
towards the substrate while the Pc skeleton moves away, see Fig. 4(b). 
The different heights of the atoms relative to the substrate lower the 
symmetry upon adsorption. With vdW interaction we obtain atomic 
shifts of 0.31 A for Zn and N, 0.43 A for Zn and C, and 0.55 A for Zn 
and H, while without vdW interaction we have 0.36 A for Zn and N, 
0.56 A for Zn and C, and 0.71 A for Zn and H. These values indicate 
that the molecular distortions are underestimated when the vdW 
interaction is included^\ An increased molecule-substrate distance 
and flattened molecule, see Fig. 4(c), correspond to an underesti- 
mated interface dipole. Without vdW interaction, our results show 
that mainly the Zn-N bond is modified by the interaction with the 
substrate. As compared to the free standing molecule, the Zn-N bond 
length increases from 2.00 A to 2.02 A and the N-Zn-N bond angle is 
heavily reduced from 180° to 163°. The fact that a distortion alters the 
interface dipole can be problematic for active organometallic inter- 
faces. To counteract the effect, fluorination of the molecule is 
expected to reduce the vertical distance and consequently the 
molecular distortion. 

To understand the electronic structure of the ZnPc/Cu(lll) sys- 
tem, we evaluate the total and partial DOSs for the most stable con- 
figurations, see Fig. 5, and compare with the free standing molecule 
(Fig. 2). Due to the interaction with the Cu substrate, the electronic 
states near Ep change significantly and the band gap disappears. 
Hence, the ZnPc molecule shows a metallic behavior, see panels 
(a) in Fig. 5. It is clear from panels (b) in Fig. 5 that the Zn states 
shift to lower energy by about 1 eV as compared to the free standing 
molecule and therefore play no longer a role near Ep. Panels (c) and 
(d) in Fig. 5 show that the C Ip^ and N Ip^ states around Ep are 
suppressed and also shift to lower energy. In the energy range from 
— 8 eV to —6 eV we mainly find hybridized Zn and C 2p states 
(/>(i-hybridization), with small contributions of the N Ipy and Ipz 




Figure 4 | Structure of ZnPc on Cu(lll). (a) before relaxation, (b) after relaxation without vdW interaction, and (c) after relaxation with vdW 
interaction. 
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Figure 5 | (a) ZnPc, (b) Zn, (c) C, (d) N, (e) H, and (f) Cu partial DOSs (per atom) of ZnPc/Cu(l 1 1) obtained by the (left) GGA without vdW and (right) 
GGA with vdW approaches. 



states, which also shift to lower energy as compared to the free 
standing molecule. 

The most significant changes in the charge distribution after 
adsorption can be described as charge transfer from the Cu surface 
to the region between molecule and substrate, which comes along 
with a significant modification of the dipole perpendicular to the Cu 
surface. As a consequence of additional charge transfer within the 
ZnPc molecule, the LUMO becomes partially occupied, giving rise to 
a metallic nature. A similar behavior has been reported by Amsalem 
and coworkers^ for ZnPc on Ag(llO), which we explain by the fact 



that the work function of Ag (4.73 eV) is similar to that of Cu 
(4.7 eV). A metallic state is also in agreement with previous calcula- 
tions for pentacene^^ and benzene^^ on a Cu substrate. In general, 
partial occupation of orbitals in MPc systems can lead to extraord- 
inary electronic properties, including superconductivity^. Such mod- 
ifications are not possible by simple physisorption of the molecule on 
the substrate. 

Figure 5 shows that Cu dominates the DOS in the energy range 
from —5.2 eV to Ep, compare panles (c), (d), and (f). The peaks 
below —1.8 eV seen in panels (a) also appear in the HOMO of the 




Figure 6 | Charge density difference isosurfaces of ZnPc/Cu( 111) obtained by the (a) GGA without vdW, (b) GGA with vdW, and (c) PBEO approaches. 

Red and blue colors represent charge depletion and accumulation, respectively. The isovalue is set to ±0.005 A"^ 
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photoemission spectra of Ref. 21. On the other hand, the small peaks 
near Ep are Shockley surface states, in agreement with the experi- 
mental situation. In general, interface states determine the charge 
transfer and give rise to a complex electronic behavior. In order to 
analyze the effects on the electronic properties, we compare results 
without and including the vdW interaction in the two panels of Fig. 5. 
Switching on the vdW interaction hardly modifies the electronic 
states. Only the C and N 2py interface states are slightly shifted to 
higher energy, while the Cu DOS stays essentially the same. 

To visualize the charge transfer, we show in Fig. 6 the charge 
density difference Ap = PznPc/cu(iii) " PznPc " Pcu(iii). where 
PznPc/cu(iii)» PznPc ^nd Pcu(iii) ^re the charge densities of the full 
system, the isolated ZnPc molecule, and the Cu substrate, respect- 
ively. In Fig. 6 results are given for GGA calculations without and 
with vdW interaction as well as for the PBEO approach. The regions 
of charge accumulation (isovalue + 0.003) are displayed in blue color, 
the regions of electron depletion (isovalue —0.003) in red color. We 
observe a significant variation of the charge density at the interface 
due to the adsorption of the molecule. A strong charge accumulation 
is found just below the center of the molecule, while regions of 
depletion appear both on the lower side of the molecule (facing the 
substrate) and on the surface of the Cu substrate. The charge accu- 
mulation at the interface is a consequence of the Pauli repulsion^^'^^. 
A quantitative analysis demonstrates that about 0.09 electrons are 
transfered into the interstitial region, while about 0.03 electrons are 
displaced within the ZnPc molecule, essentially from the outer 
regions towards the Zn atom. This redistribution of charge is 
responsible for the formation of a dipole pointing from the substrate 
to the molecule. Under inclusion of the vdW interaction the calcu- 
lated charge transfer is almost the same as without vdW interaction, 
compare panels (a) and (b) of Fig. 6. Only the magnitude is slightly 
smaller, due to the overestimation of the distance between substrate 
and molecule. Although the vdW interaction affects this distance, we 
find little influence on the charge transfer and, hence, on the elec- 
tronic structure. Results from PBEO hybrid calculations are shown in 
panel (c) of Fig. 6. While we observe slightly more charge redistribu- 
tion within the substrate, as compared to the GGA results in panels 
(a) and (b), the pictures are similar otherwise. We therefore can 
conclude that in the present system the GGA gives reliable results 
for the charge transfer. 

Discussion 

Density functional theory has been used to investigate the adsorption 
geometry and electronic properties of a ZnPc molecule on Cu(lll). It 
is found that the most stable adsorption configuration is character- 
ized by minimization of the Zn-Cu distance. Our structural optimi- 
zations leave no doubt that the Zn atom is pulled to the substrate, 
while the Pc skeleton moves away. We have demonstrated that inclu- 
sion of the vdW interaction does not improve the relaxation results 
for the present system. Modification of the chemical bonding 
(mainly the Zn-N bonds) alters the energy levels of the ZnPc mole- 
cule. It turns out that charge transfer within the ZnPc molecule 
partially fills the LUMO and induces metallic states. Most impor- 
tantly, a shift of charge from the Cu(lll) substrate and the ZnPc 
molecule into the region between them is observed. Since this will 
have strong effects on the interface dipole, the charge transfer gov- 
erns the properties of this organometallic interface rather than the 
structural effects proposed in Ref 21. Hybrid calculations result in a 
similar amount of charge transfer as obtained within the GGA. 

Methods 

The calculations are performed in the framework of density functional theory and the 
generalized gradient approximation (GGA) using the Vienna Ab-initio Simulation 
Package^^. Exchange correlation effects are treated by the Perdew- Wang functionaP^. 
For comparison, the hybrid functional approach in the PBEO flavor is employed. In 
order to determine the electronic properties of a free standing ZnPc molecule, a 
geometrical optimization is performed in a 24 A X 24 A X 7 A supercell, where the 



molecule is located in the xy-plane. Both the free standing molecule and the ZnPc/ 
Cu( 111) system are treated by a full structural optimization until all the residual forces 
on the atoms are less than 0.02 eV/A. We report results obtained both without and 
under inclusion of the van der Waals (vdW) interaction. In the ZnPc/Cu(l 1 1) 
supercell, the molecule sits on a Cu( 111) substrate of 9 A thickness and is followed by 
a 16 A thick vacuum layer to avoid artificial interaction in the slab model. Our final 
ZnPc/Cu(lll) system comprises 543 atoms (57 of ZnPc and 486 of Cu). Gamma 
point sampling is employed and the electronic wave function is expanded in plane 
waves with a kinetic energy cutoff of 400 eV. Moreover, a Gaussian smearing with a 
width of 0.05 eV is implemented. 
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